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Abstract 

A  combined  creep  and  strength  model  has  been  developed  for  the  entire 
(primary,  secondary  and  tertiary)  creep  and  the  long-term  strength  of 
frozen  soil  under  multiaxial  sfress  at  both  constant  stresses  and  constant 
strain  rates  by  a  single  (unified)  consfifutive  equation.  Secondary  creep 
is  assumed  to  be  an  inflection  point  of  a  creep  curve  defining  fime  to 
failure.  Secondary  creep  rate  is  described  by  a  new  flow  law,  fhe  stress 
function  of  which  includes  fhe  first  invariant  of  the  stress  tensor.  The 
model  consists  of  four  principal  elements:  a  constitutive  equation,  a 
viscous  flow  equation  and  a  yield  criterion,  all  united  by  a  time-to-failure 
function.  The  yield  criterion  is  selected  either  in  the  form  of  a  parabolic 
(extended)  von  Mises-Drucker-Prager  model  or  a  parabolic  (extended) 
Mohr-Coulomb  rupture  model.  The  criteria  take  into  account  that,  at  a 
certain  magnitude  of  the  mean  normal  stress  (a^Qx)'  strength 

of  frozen  soil  reaches  a  maximum.  The  yield  criteria  are  included  in  the 
time-to-failure  function,  the  shape  parameters  of  which  are  independent 
of  the  loading  regime.  The  model  has  been  verified  using  test  data  on 
creep  and  the  long-term  strength  of  frozen  soil  under  triaxial  compression 
at-10°C. 


For  conversion  of  SI  units  to  non-SI  units  of  measurement  consult  ASTM  Standard 
E380-93a,  Standard  Practice  for  Use  of  the  International  System  of  Units, 
published  by  the  American  Society  for  Testing  and  Materials,  1 91 6  Race  St., 
Philadelphia,  Pa.  19103. 
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was  provided  by  the  Office  of  the  Chief  of  Engineers  under  DA  Project 
4A762784AT42,  Design,  Construction  and  Operations  Technology  for  Cold  Regions, 
Task  BS,  Work  Unit  D02,  Investigation  and  Design  for  Foundations  in  Cold  Regions, 
Technical  review  was  provided  by  F.D.  Haynes  (CRREL)  and  Dr.  R.M.W. 
Frederking  (Cold  Regions  Science  and  Technology).  The  author  expresses  spe¬ 
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comments. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional 
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Creep  and  Strength  of 
Frozen  Soil  Under 
Triaxial  Compression 

ANATOLY  M.  FISH 


INTRODUCTION 

It  is  well-known  that  the  mechanical  proper¬ 
ties  of  frozen  soil  are  similar  to  those  of  ice. 
Therefore,  constitutive  equations  and  strength 
criteria  developed  for  one  material  are  used  in¬ 
terchangeably  for  the  other.  A  considerable  con¬ 
tribution  in  developing  constitutive  laws  and  fail¬ 
ure  criteria  of  frozen  soils  and  ice  under  a  complex 
stress-strain  state  was  made  by  a  number  of 
researchers  (Vyalov  et  al.  1963,  Sayles  1973, 
Zaretsky  and  Vyalov  1971,  Klein  and  Jessberger 

1979,  Vyalov  and  Slepak  1988,  Domaschuk  et  al. 
1991,  Zaretsky  1993,  Puswewala  and  Rajapakse 
1993,  and  many  others).  A  detailed  consideration 
of  the  constitutive  models  of  frozen  soils  and  ice 
can  be  found  in  Puswewala  and  Rajapakse  (1987). 
The  basic  assumptions  on  which  these  models 
are  based  can  be  divided  into  two  groups.  The 
models  of  the  first  group  are  based  upon  the 
following  assumptions: 

1.  Creep  and  yield  (strength)  of  frozen  soil  are 
two  separate,  loosely  related  phenomena  de¬ 
scribed  by  different  constitutive  equations  and 
failure  criteria. 

2.  Creep  of  frozen  soil  is  a  three-stage  process, 
comprising  the  primary,  secondary  and  tertiary 
creep  stages. 

3.  Failure  of  frozen  soil  takes  place  at  the  end 
of  the  secondary  or  the  tertiary  creep  stages. 

A  three-dimensional  combined  creep  and  yield 
model  of  frozen  soil  presented  below  was  origi¬ 
nally  developed  by  the  author  (Fish  1991b,  1992, 
1993)  to  describe  the  rheological  behavior  of 
polycrystalline  ice  imder  multiaxial  stress.  The 
model  is  a  generalization  for  a  complex  stress- 
strain  state  of  a  unified  creep  and  strength 
model  of  frozen  soil  under  uniaxial  compression, 
developed  earlier  by  the  author  (Fish 

1980,  1983),  based  upon  substantially  different 
premises: 

1,  Creep  and  failure  of  frozen  soils  are  a  uni¬ 
fied  process  that  can  be  described  by  a  single 
constitutive  equation  that  includes  a  strength  cri¬ 
terion  of  the  material. 


2.  Short-term  creep  of  frozen  soil  can  be  con¬ 
sidered  as  a  two-stage  process,  with  the  primary 
(decelerating)  and  the  tertiary  (accelerating)  creep 
stages.  Secondary  creep  is  not  considered  a  stage 
but  only  a  (transition)  point  (M)  on  the  creep  curve 
(Fig.  1)  defining  time  to  failure. 

3.  The  process  of  frozen  soil  failure  begins  at 
the  start  of  deformation  immediately  after  ap¬ 
plication  of  the  load  and  takes  place  at  all  stages 
of  creep. 

Studies  have  shown  that  constitutive  laws  and 
failure  criteria  based  upon  the  above  premises 
led  to  more  rigorous,  complete  and  accurate 
constitutive  laws  and  failure  criteria  for  frozen 
soil  and  ice  and  to  a  successful  solution  of  a  variety 
of  practical  engineering  problems  in  cold  regions 
(Puswewala  and  Rajapal^e  1993,  and  others).  The 
present  model  consists  of  four  principal  elements: 
a  constitutive  equation,  a  viscous  flow  equation 
and  a  yield  criterion,  all  united  by  a  time-to-failure 
function.  In  the  framework  of  this  model,  time 
to  failure  performs  a  crucial  role:  it  not  only 
unites  all  stages  of  creep,  which  makes  it  pos¬ 
sible  to  describe  the  entire  creep  process  under 
both  constant  stresses  and  constant  strain  rates 
by  means  of  a  single  constitutive  equation,  but 
it  also  combines  the  long-term  strength  and  the 
yield  criteria  so  that  the  latter  can  be  incorpo¬ 
rated  into  the  constitutive  equation  of  frozen  soil. 

CREEP  MODEL 

Constitutive  equation  (tj  =  const.) 

The  entire  deformation  process,  which  includes 
primary,  secondary  and  tertiary  short-term  creep 
of  frozen  soil  under  multiaxial  stress  at  constant 
temperature,  can  be  presented  as  a  product  of 
the  flow  equation  and  the  nondimensional  time 
function  (Fish  1980, 1983, 1991b,  1993;  Fish  and 
Assur  1984) 

where 


Triaxial 

Creep 

Tests 


Triaxial 

Strength 

Tests 


t  =  t/t^  =  Y/Y|m»  Normalized  Time  (shear  strain) 


Figure  1.  Creep  model  (after  Fish  1983,  1991b), 


yft)  =  2}l^^=  shear  strain  rate  (intensity) 
Y^=  minimum  shear  strain  rate 


/2e-  ”[{^1  "£2)^  +(^2  +(^3 


(2) 


where 

he  ~  second  invariant  of  the  deviatoric  strain 
rate  tensor 

ei  2  3  =  principal  strain  rates. 


The  nondimensional  time  function  is 

( f-iY 

F{t)  =  \  ^  =  exp{5(f  -  In i  - 1)} 


where 


(3) 

(4) 


t  =  normalized  time 
=  time  to  failure 

5  =  dimensionless  parameter;  0  <  8  <  1. 


The  function  F(t)  is  such  that  for  primary  and 
tertiary  creep,  f(f )  >  1,  and  for  secondary  creep 


(points  M  in  Fig.  1),  f  =  F  =  1  and  F(  t) 

=  1. 

The  parameter  5  is  the  first  (shape)  param¬ 
eter  of  a  creep  curve.  For  frozen  soil  5  is  between 
0.3  and  0.9,  with  an  average  value  of  --  0.6  -  0.7. 
The  effect  of  changes  of  5  on  the  shape  of  creep 
curves  is  considered  in  Fish  (1987,  1989). 

It  is  essential  that  the  creep  process  be  con¬ 
sidered  not  in  terms  of  real  time  but  in  terms 
of  normalized  time,  F  =  t/t^,  where  is  the  time 
to  failure  (durability)  of  frozen  soil,  dependent 
upon  the  type  of  soil,  its  structure,  the  unfro¬ 
zen  water  (ice)  content,  the  temperature,  the 
applied  stress,  the  loading  regime  and  other 
factors.  Then,  all  the  shape  variations  of  the  creep 
curves  caused  by  the  stress  or  temperature 
changes,  or  both,  are  expressed  by  variations  of 
a  single  parameter,  8.  In  the  certain  stress  do¬ 
main,  parameter  8  may  be  assumed  to  be  inde¬ 
pendent  of  stress  (and  temperature).  Then,  families 
of  creep  (and  stress-strain)  curves  generated  over 
a  wide  spectrum  of  mean  normal  stresses  and 
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(7) 


various  temperatures  can  be  superimposed  on 
a  single  plot  similar  to  that  in  Figure  1.  One  of 
the  principal  purposes  of  the  author  in  devel¬ 
oping  this  creep  model  (Fish  1980)  was  to  de¬ 
scribe  the  entire  creep  process  at  a  constant  stress 
and  a  constant  strain  rate  under  a  simple,  as  well 
as  a  complex,  stress-strain  state  by  means  of  one 
parameter  (8)  only. 


4>(f)  =  am/for^^"  ; 

where 

F=  ~  nondimensional  time  to  failure 
=  relaxation  time  (see  eq  12  below) 
n  =  dimensionless  parameter;  n  >  1. 

By  combining  eq  6  and  7,  the  failure  or  the  creep 
strength  criterion  of  frozen  soil  takes  the  form 


Flow  equation 

The  minimum  shear  strain  rate  is  related  to 
the  time  to  failure  by  means  of  the  flow  equa¬ 
tion  (Fish  1980) 

Yim  ”7  Yim^m 

where  C  is  the  second  shape  parameter  of  a  creep 
curve,  which  is  equal  to  the  viscous  strain  at 
failure.  Both  5  and  C  are  assumed  to  be  constant 
for  a  given  soil  temperature.  The  tilde  (--)  indi¬ 
cates  that  in  the  general  case  C  may  be  a 
function  of  stress  and  temperature  (Fish  1987, 
1989).  The  magnitudes  of  C  are  between  0.01  and 
0.3.  A  confirmation  of  validity  of  eq  5  for 
test  data  on  the  uniaxial  compression  of  ice 
and  frozen  soil  can  be  found  in  Mellor  and  Cole 
(1982, 1983),  Fish  (1983),  Gonze  et  al.  (1985)  and 
others. 

Note  that  the  mathematical  structure  of  eq  5, 
which  is  fundamental  for  this  model,  is  reminis¬ 
cent  of  the  equation  used  by  Garofalo  (1965)  and 
many  others.  However,  eq  5  is  substantially 
different  from  the  latter.  In  Garofalo's  equation, 
the  product  is  assumed  to  be  a  constant  value 
of  the  minimum  strain  rate  by  the  time  to  rup¬ 
ture,  (where  is  the  time  of  total  soil  failure 
at  the  end  of  the  tertiary  creep  stage). 

Time  to  failure 

When  time  t  =  t^\he  strain  rate  reaches  a  mini¬ 
mum,  Yi=  Yiin(Yi  =0)/  it  is  assumed  that 
failure  occurs.  The  magnitude  of  shear  stress  that 
causes  failure  at  time  can  be  expressed  (Fish 
1991b)  as  a  product  of  two  independent  func¬ 
tions — the  yield  function  (criterion)  and  the 
nondimensional  time  function  0(?) 

l2=J2o^Hi)  (6) 

where 

/2o  =  "instantaneous"  yield  criterion  (Fig,  2) 

J2  =  second  invariant  of  the  deviatoric  stress 
tensor. 

The  time  function  can  be  selected  in  the  simplest 
form  (Fish  1983): 


or 


n 


2 


(8) 

(9) 


where  x- =  and  can  be  defined 

as  the  shear  stress  and  the  "instantaneous"  shear 
stress  intensity  (resultant),  respectively  (see  eq 
23  and  26  below). 


New  flow  law 

Combining  eq  5, 8  and  9  gives  the  flow  equation 
of  frozen  soil  in  a  multiaxial  stress  state  in  the 
form 


or 


n 


V^io 


(10) 


Yim 


-E/RT 

■  h 


in  which 


(11) 


(12) 


~  FrenkeFs  relaxation  time  (Frenkel  1947) 
h  =  Planck's  constant 
k  =  Boltzmann's  constant 
T  =  absolute  temperature 
E  =  activation  energy 
R  =  gas  constant. 

Both  parameters  n  and  tJE)  are  strongly  tem¬ 
perature-dependent.  Temperature  dependencies 
of  these  parameters  are  discussed  in  Fish 
(1985). 

Note  that  the  mathematical  structure  of  the 
stress  function  is  similar  to  that  of  the  popular 
Norton-Glen  power  law  used  by  many  investi¬ 
gators  (Vyalov  et  al.  1963,  Sayles  1968  and 
others).  Despite  the  fact  that  eq  10  and  11  con¬ 
tain  a  power  function  of  stress,  they  are  hmda- 
mentally  different  from  the  well-known  Norton- 
Glen  (1958)  flow  law.  In  contrast  to  the  Norton- 
Glen  equation,  eq  10  and  11  contain  param¬ 
eters  C  and  which  have  a  definite  physical 
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=  ^  =  Mean  Normal  Stress 


Figure  2.  Yield  criteria  of  frozen  soil  and  ice:  1  von  Mises;  2  Drucker-Prager  (Drucker  and  Prager  1952);  3 
Parabolic  (Fish  1991a,  1993), 


T|  *  }/J2 


Figure  3.  Creep  strength  criterion  of  frozen  soil  and  ice  (Fish  1991a,  1993), 


meaning,  and  a  denominator  of  the  stress  func-  yield  criterion  into  the  time-to-failure  function 
tion  that  is  a  temperature-dependent  yield  cri-  makes  it  possible  to  relate  the  minimum  shear 

terion.  When  the  mean  normal  stress  strain  rate  and  the  shear  stress  over  the  entire 

the  yield  criterion,  which  is  a  function  of  the  spectrum  of  confining  pressures  (mean  normal 

first  invariant  of  the  stress  tensor  (see  eq  23  stresses).  The  validity  of  eq  10  is  confirmed  (Fish 

below),  serves  as  a  depressant  for  the  strain  rate.  1992)  by  test  data  on  secondary  creep  of  poly- 

When  the  mean  stress  the  yield  func-  crystalline  ice  imder  triaxial  compression  (Jones 

tion  in  the  denominator  of  eq  10  and  11  accel-  and  Chew  1983). 

erates  the  creep  strain  rate.  Introduction  of  the  At  a  high  stress  level,  when  tj  =  tjo  and 
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=  eq  10  becomes 


Yio=- 


(13) 


where  can  be  defined  as  the  '"instan¬ 

taneous"  shear  strain  rate,  i.e.,  the  rate  at  which 
the  transition  from  the  brittle  to  the  ductile  mode 
of  failure  takes  place. 


Creep  shear  strain 

Creep  shear  strains  are  obtained  by  integrat¬ 
ing  eq  1.  The  entire  creep  process  as  described 
by  Fish  (1983)  is 

Yi(t)  =  ^  +  YitV  (14) 

where 

X.  =  shear  stress 

G  =  shear  modulus 

yj^=  =  creep  shear  strain 

+(£2-£3f +(£3-£lf]  = 

second  invariant  of  the  deviatoric  strain  tensor 
£^23  =  principal  strains 
'  \|/  =  \|/(F)  =  integration  coefficient  calculated 
depending  upon  the  normalized  time. 

Strain  at  failure 

Failure  occurs  when  t  =  t  =  1,  and  Yi  = 
Yim-  Equation  14  gives  an  expression  for  failure 
shear  strain  (Fig.  1) 


tude  of  strain  at  failure  (point  M  in  Fig.  1)  re¬ 
tains  a  constant  value.  On  the  other  hand,  test  data 
on  triaxial  compression  of  frozen  soil  (Gorodetskii 
1975)  show  that  magnitudes  of  the  failure  strain 
vary  with  the  change  of  confining  pressure. 


Creep  strain  calculations 

Calculations  of  creep  deformations  are  con¬ 
siderably  simplified  by  using  an  approximate  in¬ 
tegral  of  eq  14  (Fish  1983  [eq  38],  Fish  1987  [eq 
20]): 


r  6i(F-i) 
+  e 


(18) 


where  -  X  and  X  is  also  a  shape  param- 

eter  of  a  creep  curve.  When  t=t^,  t  =  1  and  eq 
18  becomes  eq  16. 

Thus,  the  entire  creep  process  of  frozen  soil 
in  a  multiaxial  stress  state  at  constant  shear  stresses 
is  described  by  means  of  only  two  parameters, 
C  and  5(X-),  to  be  determined  from  test  data.  An 
evaluation  procedure  for  these  parameters  is  pre¬ 
sented  in  Fish  (1983,  1987,  1989). 


Triaxial  constant  strain  rate  tests  (  Yi  =  const) 
The  strength  of  frozen  soil  is  determined  from 
tests  in  which  the  shear  strain  rate  is  maintained 
constant  (yi  =  Yim);  then  the  shear  stress  becomes 
a  function  of  this  rate  (Fish  1983).  Defining 
Yi  =  Yif.  Yim  =  Yim'm/  C  =  and  combining  the 
latter  with  eq  10  gives 


Yim=^  +  Yim^mVm  (1^) 

where  \|f is  an  integration  coefficient  calculated 
for  F  =  1.  The  first  term  in  the  right  side  of  eq 
15  is  the  instantaneous  strain  ignored  in  Figure  1. 
Combining  eq  5  and  15  gives 

(16) 

where  ~1/(1-S)^'^^  and  X  =  1-5. 

Failure  strain  criterion 

For  frozen  soil  the  instantaneous  strains  are 
small  and  in  many  instances  can  be  ignored,  sub¬ 
stantially  simplifying  the  parameter  evaluation 
procedure.  Indeed,  in  this  case  eq  16  becomes 

or  X  =  [c/yl^f  (17) 

It  is  interesting  to  note  that  eq  17  may  be  used 
as  a  strain  failure  criterion.  Creep  test  data  from 
the  imaxial  compression  of  polycrystalline  ice 
(Mellor  and  Cole  1982)  showed  that  the  magni¬ 


r  —  T 

”Uln,J 


Nl/n 


(19) 


Subscript  tn  indicates  that  the  magnitudes  of 
stress  (x^j^)  and  strain  (y^^)  are  referred  to  the 
point  M  corresponding  to  the  maximum  (peak) 
shear  stress  in  Figure  1.  In  the  particular  case 
when  the  applied  strain  rate  Yi  =  Yim  =  Yio' 

Yim  =  Yio  =  Yio'o  the  peak  shear  stress  11;^  = 
XjQ  is  defined  as  the  "instantaneous"  yield  stress 
(criterion)  of  frozen  soil  and  ice  and  is  given  by 
eq  23  (26). 

Introduction  of  the  yield  criterion  Xj^  and  pa¬ 
rameter  into  eq  19  makes  the  latter  essentially 
different  from  similar  equations  of  Nadai  (1950), 
Sayles  and  Epauchin  (1966),  Haynes  et  al.  (1975), 
Ladanyi  (1972),  Baker  (1979),  and  others.  Cor¬ 
respondence  between  constant  stress  (CS)  creep 
tests  and  constant  strain  rate  (CSR)  tests  for 
uniaxial  compression  of  ice  and  frozen  soil  is 
confirmed  by  the  analysis  of  test  data  (Mellor 
and  Cole  1983,  Fish  1983  and  others).  The  va¬ 
lidity  of  eq  19  to  test  data  on  the  strength  of 
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polycrystalline  ice  under  triaxial  compression  in 
a  wide  spectrum  of  strain  rates  is  shown  in  Fish 
(1991b,  1992,  1993). 

Stress-strain  relationships  outside  point  M  are 
obtained  by  combining  and  rearranging  eq  1,  3 
and  19 


(20) 


or 


/  ^5/n  r 

Yi 

exp 

UimJ 

Yi 


-1 


V  I  im 


(21) 


where  t  =  t/t^=yi/yi^._ 

For  secondary  creep,  F  =  1  and  the  strain  rate 
dependency  of  stress  is  defined  by  eq  19.  The 
analysis  of  CSR  test  data  showed  (Fish  1983)  that 
eq  21  gives  a  satisfactory  approximation  for  the 
stress-strain  relationships  of  frozen  Fairbanks  silt 
in  uniaxial  compression. 


YIELD  MODEL 

Studies  (Chamberlain  et  al,  1972,  Sayles  1973, 
Gorodetskii  1975,  Parameswaran  and  Jones  1981, 
Baker  et  al.  1982,  Jones  1982,  Jones  and 
Parameswaran  1983)  show  that  the  shear  strength 
of  frozen  soil  and  ice  imder  triaxial  compression 
is  a  nonlinear  function  of  confining  pressure. 
Moreover,  at  a  certain  magnitude  of  confining 
pressure  defined  (Chamberlain  et  al.  1972)  as  ice 
melting  pressure,  the  shear  strength  of  frozen  soil 
(and  ice)  reaches  a  maximum  value,  A  further 
increase  of  the  pressure  leads  to  a  decrease  of 
the  frozen  soil  (ice)  strength.  On  the  other 
hand,  in  the  small  stress  domain,  a  linear 
relationship  is  observed  between  the  shear 
strength  of  frozen  soil  and  the  applied  confin¬ 
ing  pressures. 

A  detailed  consideration  of  yield  criteria  ap¬ 
plied  to  frozen  soil  can  be  found  in  Vyalov  et 
al.  (1963, 1981),  Vyalov  (1986),  Sayles  (1973, 1988), 
Fish  (1991a)  and  others. 

In  general  terms  a  yield  (strength)  criterion 
for  frozen  soil  (ice)  imder  multiaxial  stress  at  con¬ 
stant  temperature  that  satisfies  the  above  test  data 
can  be  expressed  in  the  form 


(22) 


where 

=  first  invariant  of  the  stress  tensor 


/2  =  second  invariant  of  the  stress  deviator. 
Studies  show  that  a  parabolic  yield  criterion  (Fig. 
2)  developed  by  Fish  (1991a)  well  describes  test 
data  on  the  strength  of  frozen  soil  and  ice 
over  a  wide  range  of  confining  pressures  and 
various  applied  strain  rates  or  times  to  failure, 
i.e. 


Tj  =  c  +  ba^ 


b  2 

^^max 


where 


(23) 


+  (02  -  Oaf  +(03 
=  shear  stress  (intensity) 


a 


m 


1 

3 


(01  +  02+03)-  mean  normal 


stress 

01,02,03  =  principal  stresses 

c  =  resistance  in  pure  shear  associated  with 
the  frozen  soil  cohesion  in  the  octahe¬ 
dral  plane 

b  =  tan()),  where  ([)  =  angle  of  internal  fric¬ 
tion  in  the  octahedral  plane 
0max  =  melting  pressure,  i.e.,  the  magni¬ 
tude  of  the  mean  normal  stress  at 
which  the  shear  stress  of  frozen  soil 
reaches  a  maximum 


max  Tj  = 


b 

^  +  2  ^niax 


(24) 


The  yield  curve  intersects  the  hydrostatic  axis 
at  points  Hi  2/  the  abscissas  of  which  are 


^1,2  “  ^max  — 


2c 


lU2 


^max  ^  ^  ^max 


(25) 


Note  that  the  yield  curve  shape  is  defined  by 
the  ratio  of  parameters  c  and  b. 

Equation  23  can  be  considered  as  an  extended 
von  Mises-Drucker-Prager  yield  criterion.  At  low 
stress  levels  (Omax  ^°°)/  23  transforms  into 

the  Drucker-Prager  (1952)  yield  criterion.  For 
frictionless  materials  {b  =  0),  eq  23  reduces  to  the 
von  Mises  yield  criterion.  Parameter  ^^7 
be  regarded  as  one  of  the  fundamental  physi¬ 
cal  characteristics  of  frozen  soil,  probably  closely 
related  to  the  phase  changing  pressure  (point  Hi) 
for  a  given  temperature  (Nadreau  and  Michel 
1986).  Point  H2  with  abscissa  h2  may  be  considered 
as  the  resistance  of  frozen  soil  to  triaxial  exten¬ 
sion. 

Equation  23  describes  the  strength  of  frozen 
soils  in  a  multiaxial  stress  state,  irrespective  of 
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the  time  to  failure  or  the  applied  strain  rate  by 
three  temperature-dependent  parameters — c,  b((j)) 
and  — all  of  which  have  definite  physical 

meanings  and  are  easily  determined  from  test 
data.  It  can  be  shown  that  a  parabolic  (extended) 
Coulomb-Mohr  yield  criterion  (Fish  1991a)  can 
also  be  applied  to  describe  such  a  peculiar  be¬ 
havior  of  frozen  soil  under  a  high  confining 
pressure. 


CREEP  STRENGTH 


Creep  strength  criteria 

It  has  been  shown  above  that  the  creep  strength 
of  frozen  soil  is  a  function  of  time  or  applied  strain 
rate.  The  absolute  values  of  the  frozen  soil 
strength,  as  well  as  the  magnitudes  of  the  strength 
parameters,  change  from  the  greatest  (instanta¬ 
neous)  values  at  the  time  of  loading  (f^)  to  zero 
when  tin  Consequently,  the  creep  strength 
of  frozen  soil  is  characterized  by  a  family  of  curyes 
for  given  times  fo  <  given  strain  rates 

YiO>Yii>Yin  (Fig-  3)- 

A  grapho-analytical  procedure  for  evaluation 
of  the  frozen  soil  strength  based  upon  triaxial 
creep  and  strength  tests  data  has  been  developed 
by  Sayles  (1973)  using  the  Coulomb-Mohr  rup¬ 
ture  model. 

An  analytical  procedure  for  interpretation 
of  triaxial  creep  and  constant  strain  rate 
test  data  and  determination  of  strength  param¬ 
eters  of  frozen  soil  (and  ice)  has  been  devel¬ 
oped  by  Fish  (1991a,  b,  1993,  1994)  based 
upon  the  parabolic  (extended)  creep  strength 
criteria. 

It  has  been  shown  (Fish  1991a)  that  the  creep 
strength  of  frozen  soil  can  be  expressed  as  a 
product  of  the  parabolic  yield  criterion  (eq  23) 
and  a  nondimensional  time-to-failure  function 
(eq  7),  i.e. 


'fi(0  =  XioO(0  = 


2an 


\m 

(26) 


Equation  26  is  a  parabolic  criterion  for  creep 
strength.  When  0(t)  =  1  and  eq  26  be¬ 

comes  eq  23.  The  subscript  0  indicates  that  the 
failure  shear  stress  or  the  magnitudes  of  the 
strength  parameters  refer  to  the  instantaneous 
yield  condition.  Accordingly,  the  maximum  shear 
strength  is  defined  by 

maxTi(f)  =  ^Co  +  y Omax  •  (27) 


One  may  assume  that  the  ratio  of  parameters 
b  and  c  in  eq  25  is  constant,  i.e.,  the  positions 
of  points  and  H2  are  are  unaffected  by  the 
change  of  the  time  to  failure  or  the  applied  strain 
rate. 

In  the  particular  case  in  the  range  of  small 
stresses  when  0,  the  parabolic  criterion 

eq  23  becomes  the  modified  Drucker-Prager  cri¬ 
terion  for  creep  strength  (Fish  1991,  1994) 

^iit)  =  ico  +  bo<Jmy^0)  ■  (28) 

When  b  =  0,  and  eq  28  transforms  into  the 
modified  von  Mises  yield  criterion 

Ti(t)  =  Co<I>(0  .  (29) 

Equations  26  through  29  are  valid  for  constant 
strain  rate  (CSR)  tests  in  which  0(0  =  <I)(Yi),  given 
by  eq  31  and  32. 

Creep  strength  criteria  for  frozen  soil  can  also 
be  selected  in  d\e  form  of  a  parabolic  (extended) 
Coulomb-Mohr  criteria,  discussed  in  great  de¬ 
tail  in  Fish  (1991a). 

Time-dependent  strength 

Passing  vertical  planes  through  point  M  in  Fig¬ 
ure  2,  one  obtains  the  relationship  between  the 
creep  strength  of  frozen  soil  and  the  time  to 
failure  t^^.  Such  a  dependency  relates  the  instan¬ 
taneous  (yield)  strength  of  frozen  soil  determined 
from  creep  tests  at  time  (or  from  constant  strain 
rate  tests  at  the  shear  strain  rate  y^o )  with  a  creep 
strength  x^{t)  =  'ii{t^)  determined  at  times  = 
...  etc.  During  creep  at  constant  stresses, 
the  shear  strength  of  frozen  soil  as  well  as  the 
strength  parameters  in  eq  23  become  functions 
of  time  (Fig.  3),  i.e. 

=  (30) 

in  which 

Ti(0  =  XioO(0,  c(0  =  CoO(?)  and  b{t)^b^^{t) 
Comparing  eq  7,  9,  26  and  30  gives 


Xi(0_i?(t)^c(0^/.  yl/n 


(31) 


Consequently,  for  constant  strain  rate  tests,  time 
variations  of  the  strength  and  the  strength  pa¬ 
rameters  in  eq  26  through  31  are  defined  by  the 
time  function 
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-1/n 


<D(t)  =  (D(Yi)  = 


'ti(Yi)  KVi)  ■■  ^(ti) 


Tim 

i'i^oJ(32) 


given  by  eq  19. 

Note  that  the  two  parameters  and  n  define 
the  shapes  of  the  curves  of  the  long-term  (creep) 
strength  of  frozen  soil.  It  will  be  shown  below 
for  frozen  soils  at  six  different  types  of  loading 
regimes  (uniaxial  compression  and  tension,  pure 
shear,  and  triaxial  compression  at  three  differ¬ 
ent  mean  normal  stresses)  that,  for  all  practical 
purposes,  these  parameters  may  be  considered 
to  be  independent  of  the  loading  regime. 

Thus,  the  entire  process  of  deformation  and 
time-dependent  failure  of  frozen  soil  in  a  mul- 
tiaxial  stress  state  at  both  constant  stresses  and 
constant  strain  rates  is  described  by  seven  pa¬ 
rameters: 

Creep  5(X)  and  C 

Failure  time  and  n 
Yield  c^, 

which  are  determined  from  test  data. 


obtained  by  S.E.  Gorodetskii  (Vyalov  et  al.  1963, 
1979,  Gorodetskii  1975).  Cylindrical  samples  of 
Kellovian  sandy  silt  with  a  water  content  of  26%, 
a  density  of  1.81  g/cm^  and  a  specific  gravity 
of  2.71  were  used  in  the  tests.  The  test  procedures 
and  the  sample  preparation  technique  are  dis¬ 
cussed  in  great  detail  in  Vyalov  et  al.  (1963, 
1981). 

Creep  tests  were  carried  out  at  six  different 
loading  regimes:  uniaxial  compression  and  ten¬ 
sion,  pure  shear  (torsion)  and  triaxial  compres¬ 
sion  at  three  different  magnitudes  of  the  mean 
normal  stress — =  1.5,  3.0  and  6.0  MPa.  The 
tests  were  carried  out  at  a  temperature  of  0  -lO'C. 

The  magnitudes  of  shear  stresses  in  the  triaxial 
creep  tests  (02  =  CT3) 

Xi=(0l-03)/^/3 

were  maintained  constant  by  adjusting  the  mag¬ 
nitudes  of  the  axial  stresses  in  the  cross  section 
of  the  samples.  In  each  series  of  tests,  the  magni¬ 
tudes  of  the  mean  normal  stresses 


TEST  DATA 

Time-dependent  strength  of  frozen  soil 

Test  data  on  the  time-dependent  (creep) 
strength  of  frozen  soil  under  various  loading 
regimes  used  in  the  following  analysis  were 


Om=  (01  +  203)73 

were  also  maintained  constant. 

Test  data  on  the  long-term  strength  of  Kellovian 
silt  at  -lO^'C  under  various  loading  regimes  are 
presented  in  Figure  4.  In  these  tests  the  time  to 


tf,  Time  to  Failure  (hr) 

Figure  4.  Creep  strength  of  Kellovian  silt  under  various  loading  regimes  at 
-10*C.  Data  from  Gorodetskii  (1975)  and  Vyalov  et  al  (1963,  1981). 
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failure  was  calculated  at  the  end  of  the  second¬ 
ary  creep  stage.  Therefore,  in  the  notation  of  ff, 
the  subscript /is  used  instead  of  m.  It  has  been 
shown  by  Fish  (1991a)  that  all  the  time  depen¬ 
dencies  of  frozen  soil  strength  presented  in  Figure 
4  can  be  described  by  a  single  time-to-failure 
function  (eq  9),  two  parameters  of  which 

=  0.086  hr  and  n  =  9.64 

are  independent  of  the  loading  regime. 

In  Figure  5  the  dependencies  are  presented 
of  the  shear  strength  upon  the  mean  normal  stress 
for  different  moments  of  failure  time  f£  =  1  hr 
and  =  12  hr.  One  can  conclude  that  these  de¬ 
pendencies  can  be  described  by  the  parabolic 
strength  criterion  eq  23.  It  has  been  found  that 
the  magnitudes  of  the  strength  parameters  of 
Kellovian  silt  at  -10°C  are: 

Cq  =  1.98  MPa  =  1.66  ((t)^  =  59”) 

Gmax  =  6*38  MPa  h2  =  -1.1  MPa. 

A  detailed  procedure  for  evaluation  of  these 
parameters  is  given  in  Fish  (1991). 

Figures  4  and  5  contain  both  test  data  and  pre¬ 


dicted  dependencies  of  the  creep  strength  of 
Kellovian  sandy  silt  under  various  loading  re¬ 
gimes  and  various  moments  of  time.  One  can 
conclude  that  yield  and  creep  strength  criteria, 
eq  9  and  23-31,  describe  the  test  data  well. 

In  Figure  6,  time  variations  of  the  creep  strength 
as  well  as  the  strength  parameters  of  frozen  soil 
are  presented  in  the  form  of  a  dimensionless  time- 
to-failure  function  calculated  by  eq  7.  One  can 
see  that  test  points  correlate  well  with  the  rela¬ 
tionship  predicted  by  eq  31. 

From  the  above  analysis,  we  can  conclude  that 
the  shape  parameters  of  the  creep  strength  curves 
of  frozen  soil  are  independent  of  the  loading 
regime,  i.e.,  the  magnitudes  of  the  shape  param¬ 
eters  and  n  do  not  vary  with  the  change  of 
the  loading  regime. 

Creep  of  frozen  soil  under  triaxial  compression 

In  Figure  7  relationships  are  presented  both 
for  tests  and  for  predicted  creep  curves  of  fro¬ 
zen  Kellovian  silt  under  triaxial  compression 
(G2  =  <73)  at  temperature  -10°C.  Test  creep  curves 
were  obtained  by  S.E.  Gorodetskii  (Vyalov  et  al. 
1979,  Gorodetskii  1975)  for  various  shear  stresses: 

=  2.6,  2.29  and  1.73  MPa.  The  predicted  rela- 


Figure  5.  Shear  strength  vs  G^,  mean  stress  for  Kellovian  silt  at  -10"C.  Data  from 
Gorodetskii  (1975), 
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Figure  6.  Normalized  curve  of  the  long-term  strength  of  Kellovian  silt  at 
-WC.  Data  from  Gorodetskii  (1975),  Vyalov  et  al  (1963,  1981). 


tionships  were  calculated  by  eq  17  and  18  us¬ 
ing  magnitudes  of  creep  parameters:  X  =  0,2  and 
C  =  0.16, 

The  magnitudes  of  the  times  to  failure  t^  nec¬ 
essary  for  these  predictions  were  calculated  by 
eq  9  and  by  using  the  strength  parameters  of 
Kellovian  silt  at  -lO^C  obtained  above. 

One  can  conclude  that  calculated  dependen¬ 
cies  of  shear  strain  as  a  function  of  time  corre¬ 
spond  to  test  data  on  the  creep  of  frozen  soil  under 


triaxial  compression  and  can  be  used  for  solv¬ 
ing  practical  problems  in  cold  regions. 

CONCLUSIONS 

1.  The  combined  strength  and  creep  model  de¬ 
veloped  above  describes  the  entire  (primary,  sec¬ 
ondary  and  tertiary)  creep  and  the  long-term 
strength  of  frozen  soil  under  triaxial  compres- 


Figure  7.  Creep  curves  of  Kellovian  silt  under  triaxial  compression 
(02  =  03)  at  Data  from  Vyalov  et  al.  (1979). 
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sion  at  both  constant  stresses  and  constant  strain 
rates  by  a  single  (unified)  constitutive  equation. 

2.  The  three-dimensional  combined  constitutive 
model  consists  of  four  principal  elements:  a 
constitutive  equation,  a  viscous  flow  equation 
and  a  yield  criterion,  all  united  by  a  time-to-failure 
function. 

3.  The  time-to-failure  function  or  the  equation 
of  the  long-term  (creep)  strength  of  frozen  soil 
includes  a  parabolic  criterion  for  creep  strength, 
which  is  selected  in  the  form  of  the  modified 
(extended)  von  Mises-Drucker-Prager  criterion 
or  the  modified  (extended)  Coulomb-Mohr  yield 
criterion  (Fish  1991a). 

4.  Time  variations  of  the  frozen  soil  strength 
as  well  as  the  strength  parameters  are  described 
by  the  curve  of  the  long-term  (creep)  strength, 
the  shape  parameters  of  which  are  independent 
of  the  loading  regime  and  can  be  determined  from 
tests  under  a  sample  stress-strain  state  such  as 
uniaxial  compression  or  pure  shear. 

5.  The  model  has  been  verified  using  test  data 
on  the  creep  and  time-dependent  failure  of 
Kellovian  sandy  silt  at  ~10°C. 
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